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, The random forced Navier-Stokes equation can be obtained as a variational problem of a proper 

action. In virtue of incompressibility, the integration over transverse components of the fields allows 
^"*^ ' to cast the action in the form of a large deviation functional. Since the hydrodynamic operator is 

nonlinear, the functional integral yielding the statistics of fluctuations can be practically computed 
by linearizing around a physical solution of the hydrodynamic equation. We show that this procedure 
yields the dimensional scaling predicted by K41 theory at the lowest perturbative order, where the 
1 1 . perturbation parameter is the inverse Reynolds number. This result is valid over a finite spatio- 

Sh ' temporal domain, where the physical solution can be considered as stationary. 

t3 : 

1 I. INTRODUCTION 

^ , A field theoretic approach to the study of the random stirred Navier-Stokes equation (NSE) can be traced back to the 
O ' seminal paper by Martin, Siggia and Rose [1] . This was the starting point for the application of many perturbative 
! strategies, e.g. diagramatic expansions and renormalization group methods [2]. The many technical difficulties 
£>v encountered in developing such approaches avoided to gather conclusive achievements. For recent developments 
' and applications the reader can be addressed to [3]. In this paper we show that one step forward along this field- 
theoretic approach allows one to cast the action associated with the NSE into the form of a large deviation functional. 
Recently, large-deviation theory has scored sensible success in describing fluctuations in stationary non-equilibrium 
t-H , regimes of various microscopic models [4]. This approach is mainly based on the extension of the time-reversal 
conjugacy property introduced by Onsager and Machlup [5] to stationary non-equilibrium states. In practice, thermal 
0^ ' fluctuations in irreversible stationary processes can be traced back to a proper hydrodynamic description derived from 
the microscopic evolution rules. The action functional has the quadratic form 

W)]0°) = 5 r dt(W,K(p)W) (1) 

o . 

where p(t, x) represents in general a vector of thermodynamic variables depending on time t and space variables x. 
The symbol (• , •) denotes the integration over space variables. IF is a hydrodynamic evolution operator acting on 
p: it vanishes when p is equal to the stationary solution p, which is assumed to be unique. The positive kernel K(p) 
represents the stochasticity of the system at macroscopic level. According to the large deviation-theory the entropy 
C"| [ S of a stationary non-equilibrium state is related to the action functional / as follows: 



Oh- 



- 1—^ 

X 



S(p)=inf/ hoo , ](/3) (2) 



p 



where the minimum is taken over all trajectories connecting p to p. 



For our purposes it is enough to consider that the action functional / provides a natural measure for statistical 
fluctuations in non-equilibrium stationary states, so that, formally, any statistical inference can be obtained from /. 
Indeed, from the very beginning we have to deal with a hydrodynamic formulation, namely the random stirred NSE: 
in the next Section we will argue that an action functional of the form (1) can be obtained by field-theoretic analytic 
calculations. 

In particular, explicit integration over all longitudinal components of the velocity field and over the associated 
auxiliary fields can be performed. This allows to obtain a hydrodynamic evolution operator W which depends only 
on the transverse components of the velocity field vx{t,x). Moreover, the positive kernel K amounts to the inverse 
correlation function of the stochastic source, while any dependence on the form of the pressure tensor and of the noise 
does not enter in the determination of IF. Accordingly, many of the technical difficulties characterizing standard 
perturbative methods and diagramatic expansions have been removed. 
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On the other hand, we have to face with new difficulties. The hydrodynamic operator appearing in the large 
deviation functional is nonlinear, so that functional integration is unfeasible. One should identify a stationary solution 
vx{t,x) of the associated hydrodynamic equation and linearize the hydrodynamic operator around such a solution. 
Then, functional integration could be performed explicitly on the "fluctuation" field. In order to be well defined 
this approximate procedure would demand the uniqueness of the stationary solution of the nonlinear hydrodynamic 
equation. Conversely, it can be easily verified that it admits several solutions. This notwithstanding, we have found 
only one solution that does not yield unphysical divergences in the long time and large space limits (see Section III). 
Accordingly, we have assumed that fluctuations can be meaningfully estimated only with respect to this solution. 
Specifically, our statistical non-equilibrium measure is constructed by considering "trajectories" of the transverse 
component of the velocity field vr(t,x) that are connected to Vr(t,x) by the fluctuation field UT(t,x). Let us point 
out that, consistently with the linearization procedure, the statistical measure is assumed to be concentrated around 
VT(t,x). One further technical problem is that the linearized hydrodynamic operator has coefficients depending on 
both space and time variables through VT(t, x), as a consequence of its originary nonlinear nature. It depends also on 
the Reynolds number 1Z in such a way that the action integral can be naturally solved by a perturbative expansion in 
IZ^ 1 . We want to remark that sufficiently large values of 1Z also yield a weak dependence on time of VT{t,x). These 
points are discussed in Section IV. 

Since our main purpose here is the estimation of the structure function (see Section V) as an average over the non- 
equilibrium measure induced by the action /, we have to assume that the perturbative expansion applies in a wide 
range of values of 1Z. In particular, we guess that it holds also for moderately large 1Z, since a statistical average of 
any observable cannot be valid for too large values of 7Z, i.e. in a regime of fully developed turbulence. We will argue 
that statistical estimations can be consistently obtained for values of 1Z which extend up to the region of stability of 
the solution v T (t, x). Beyond this region wc have no practical way of controlling the convergence of the linearization 
procedure. It is worth stressing that we obtain an analytic expression of the structure function: the so-called K41 
scaling law [6] is recovered on a spatial scale, whose nontrivial dependence on 1Z is explicitly indicated. 



II. THE MODEL 



We consider the Navier-Stokes equation for the velocity vector-field components v a describing a divergence-free 
homogeneous isotropic flow: 



d_ 

dt 



^V 2 ) v -{x,t)+v^{x,t)^v a {x,t) + ~P(x,t) - r(x,t) = 0, 



— v a {x,t) =0. 
dx a x ' 



(3) 
(4) 



Here, the field f a represents a source/sink of momentum necessary to maintain velocity fluctuations. Customarily 
[7] , we assume f a to be a white-in-time zero-mean Gaussian random force with covariance 



(f a (x, t)f(x', t')) = F af3 (x - x') S (t - t') 



A standard choice for F is 



F{x) 



D a L 3 f 3 

= WrJ dp 



{Lpfe 



(5) 



(6) 



where D is the power dissipated by the unitary mass, L is the integral scale and the exponent s is an integer number 
of order one (a typical value is s = 2). Due to constraint (4), the field v a depends only on the transverse degrees 
of freedom of f a . Without prejudice of generality wc can also assume divergence- free forcing yielding the additional 
relation 



^(*-fO = ^(*-*) = 0. 



(7) 



Following the Martin-Siggia-Rose formalism [1] we introduce the Navier-Stokes density of Lagrangian 

-p(x,t)-r(x,t) 



L(v,w,P,QJ)=w a {x,t) 
1 d 



+ 



p dx° 



(8) 
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where the held w a is the conjugate variable to the velocity held v a and the held Q is the Lagrangian multiplier related 
to constraint (4) . In a similar way the pressure held P acts as a Lagrangian multiplier of the solenoidal constraint for 
the auxiliary held w a . The generating functional is given by the integral 

W(P,Q,J) = J VvVwVfexp^i J dt d 3 x [£(v,w,P,Q, f) + J a v a ] 

-\ J dtd^yrF-'f} 



(9) 



where J a is an external source. By performing successive integrations over the statistical measure, T>fe~ f ? F $ j on 
the auxiliary held Q, on the longitudinal components of w a and on the pressure held, all the longitudinal degrees of 
freedom can be eliminated and we end up with the effective functional for the transverse components (here denoted 
by v Ta and w Ta ) 



W(J)= / Vv T Vw T exp{i / dt d 3 x 



+W^Vj,dpVTa + J a VT 



{t,x)~ J , 



9 _2 l 

dt~ vV ]vTa 



dt d 6 xd 6 yw%(t,x)F a p(\x - y\)w£(t,y) 



(10) 



Diagramatic strategies are usually applied at this level. We want to point out that a completely different point of 
view can be followed by observing that also the transverse components of the auxiliary held u>£ can be integrated out 
and one finally obtains: 



W(J) = J Vv T e^^ I( - VT ^ +i f dtd3xJ « v T 

where the action functional I is given by 

I{v T )= j dtd 3 xd 3 y 
1 



(11) 



F^(\x-y\) 



^ - ;'V- j (■'/.(/. r) + r" r (l.r)O i r'-, (/. 7) 
d 



dt 



- vV z v p T {t, y) + i4(t, y)d x v p T {t, y) 



(12) 



This expression links the functional representation of the NSE to the large deviation theory developed in [8,4]). In 
particular the entropy is related to the functional I(vt) by [4] 



S(vt) = -inf 7(u) 
2 v 



(13) 



Where the minimum is taken over all trajectories connecting a steady-state at time t = — oo with the velocity held 
«?(0). 



III. A QUASI-STEADY SOLUTION AND ITS STABILITY 



We consider the equation of the extremal condition for the functional W( J) at J a = 0. In terms of the functional 
I(vt) it reads: 



SI{v T J) = 2 , d3y 



5vj,(t, x) 



-6 t7 a v^(t,x)dp 
+v$(t,y)d x vP(t,y) 
We observe that every solution of the equation 



1 



F^(\x-y\) 
= . 



dt 



vV 2 )v^(t,y) 



(14) 
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- iA7 2 ^ v^(t,x) +v£(t,x)d\v%,{t,x) = d p <$>{t,x) , (15) 

is solution of (14) as well. This is because of the presence of the projector on the transverse degrees of freedom in 
F a/3 (\x — y\). In particular the solutions of equation (15) are forcing independent. Indeed this equation is the usual 
unforced NSE with $ = -P. We consider now the particular case where d 13 ® = 0, that corresponds to search for 
transverse solutions satisfying the relation 

df)(y$(t,x)d x v^{t,x)) =0 . (16) 

In order to investigate the statistics of fluctuations, the determination of the hydrodynamic trajectory minimizing 
the entropy functional has to be performed explicitly. Unfortunately, existence and uniqueness theorems for the 
solution of eq. (15) are not available and the criterion for chosing a suitable solution can rely only upon physical 
considerations. Actually, we have found several different solutions: among them, the only one unaffected by divergences 
in space and time is the following: 

V«(t,x) = ^jl + e-^ t sm( , 2U _ { "f}' 3 )\ , with t>0, 
V a (t,x) = ^ , with t < . 

Notice that this solution decays exponentially in time with a rate vjL 2 = t^, 1 , where td is the diffusion time scale. 
Moreover, for <C L this solution approximates a linear shear flow. This is well known to produce small scale 
instabilities for sufficiently large Reynolds number. 

Due to the invariance under Galileo transformations, solution (17) is determined up to two vector parameters; the 
velocity U a and the rotation axis b a . These constants are related to the energy and to the enstrophy at the time t. 
Condition (16) is trivially satisfied, indeed 

V p (t,x)dpV a (t,x) = . (18) 

In other words, eq.(17) is also solution of the diffusion equation (d t — fV 2 ) u a (t,x) — 0. Obviously (17) is not a 
steady solution. It actually decays with a typical time r = 4 ffi which increases with 1Z . We assume to consider 
sufficiently large values of the Reynolds number (or, equivalently, sufficiently small values of U 2 ) so that (17) can be 
viewed as a quasi-steady solution. 

Small fluctuations around this quasi-steady solution can be analyzed by introducing the fluctuation field u a (t, x): 

v a (t,x) = V a (t,x)+u a (t,x) (19) 

and by studying the linearized equation 

^-u a {t, x) - vV 2 u a {t, x) + V p (t, x)^u a {t, x) + u?{t, x)^V a {t, x) = . (20) 

For sufficiently small initial perturbations, linear analysis can be applied up to a time scale where the nonlinear terms 
are kept small with respect to the linear ones. An upper bound for this time scale is computed in Appendix A: 

0<*<^, (2D 

This upper bound turns out to be consistently smaller than td- Moreover, a sufficient condition for the stability of 
solution (17) with respect to small perturbations has been also derived in Appendix A: 

^ > 1 ■ (22. 
This condition indicates that for high Reynolds numbers only large wave-numbers are stable. 
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IV. PERTURB ATIVE ANALYSIS OF THE GENERATING FUNCTIONAL 



By exploiting the translational invariance of the functional measure, (11) can be rewritten in the form 

W(J) = / Vu T e-htvW+if dtd 3 xJo,u% _ 



(23) 

We recall that we aim at obtaining explicit expresions for the structure functions of the perturbation field around the 
quasi-steady solution by performing derivatives of the functional generator with respect to the currents J a . 
Without any further approximation, such a program seems to be prohibitive. Some simplifications have to be intro- 
duced, yielding a structure of the functional integral which involves Gaussian integrations. 
As a first step in this direction we replace the original action in the functional (23) by the bilinear action 

I v (u)= J dtd 3 xd 3 y [(d t - vV 2 M(x) + V{x)d p v$.(x) + u p T (x)d p V a (x)] p aP{ ^_^ ) x 

\d t - W 2 y )4{y) + V\y)d x 4{y) + u*{y)d x V {y)\ + 0(u 3 ) . (24) 

which comes from the linearized equations (20) and where we have used the notation x = (t, x). Note that the extreme 
solution V a , around which fluctuations are computed, appears explicitly in the functional iy(w). We have already 
observed that the solution V a is not unique. Nevertheless, if the stability conditions of this solution with respect to 
small perturbations are fulfilled we can conclude that it cannot be influenced by other possible solutions. 

One further simplification can be introduced by observing that solution (17) naturally suggests a perturbative 
expansion in integer powers of Indeed, the Fourier transform of action (24) up to the first order in the expansion 
parameter becomes 



Iv{u) = J ^^(-^ M /(-^F^)<^ C ^ + °(^2 



(25) 



where the hydrodynamic evolution term M^ c (p)u^(p) has the form 



i \ Pa + 2? ■ u 



2 C - 

up — -p ■ U 



b A U 

—d v 

4uTZ Pt 



r , (bAXj) 
4 AvTZ 



> u c (p). 



(26) 



Here C = 



2U 



_ _ , p = (po,p) where po and p are the conjugate variables of t and x, respectively. 



^Jb 2 U 2 -(U-b) 2 ' 

In order to evaluate the functional integral in (23), the diagonalization of the matrix M p a (—p) F Jp ^ (p) is 
required. Since by definition the factor [F a @ (p)]^ 1 is proportional to the identity operator in the space of the 
transverse solutions it remains to diagonalize only the matrix M^(p). 

The computation of the eigenvalues, A, of M^(p) can be accomplished by a standard procedure, which, however 
requires lenghty and tedious calculations: they are sketched in Appendix B. We report hereafter the final form taken 
by the generating functional: 



J{H)V(t> T 



(27) 



Here H is the matrix that diagonalizes M p a (—p) F d^ M^ (p); J(H) is the Jacobian of the basis transformation 
u — ► (j), J — ► r] engendered by H. 

By performing a Gaussian integration we obtain the normalized functional in term of the r\ a source 



By returning to the representation in the original basis this equation can be rewritten as: 



W(J) 



(28) 



(29) 



The functional (29) is the starting point for the calculation of all correlation functions (and structure functions), that 
can be obtained by derivation with respect to the J£ currents. The procedure to achieve this goal is the subject of 
the next Section. 
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V. SHORT-DISTANCE BEHAVIOR OF THE SECOND ORDER STRUCTURE FUNCTION 



The expression derived for the generating functional (29), contains all the statistical information on the fluctuations 
around the basic solution V a . Here we will perform analytic calculations for the particular class of fluctuations captured 
by the lowest, nontrivial, integer moment of velocity differences between points separated by a distance r. For the 
velocity field, u a , this is the second-order structure function defined as 

52= (\u(t, f+ x) - u(t, x)\ 2 } 

= (\(u a (t,r + x)-u a (t,x))(u a (t,f+x)-u a (t,x))\), (30) 

where the brackets denote averages on the forcing statistics. 

By assuming isotropy and homogeneity of the velocity field, expression (30) is expected to assume the typical form of 
scale invariant functions 

S 2 (r)=r^F 2 (t,j) (31) 

Here r — \r\ and L is the large spatial scale associated with the noise source. It is worth stressing that, at variance with 
fully developed turbulent regimes, here the assumption of isotropy and homogeinity have to be taken as a plausible 
hypothesis allowing for analytic computations. 



We want to point out that any exponent ( n should be independent of the basis chosen for representing the functional 
W. Making use of (28), one obtains: 



S 2 (r)= (\u{x + r) -u{x)\ 2 } 



iSrj a (t, x + r) iSri a (t,x) J \i5r) a (t,x + r) iSrj a (t,x) 

dp d 3 p 



7J = 



= 2(A" a (0,f)-A" a (6))=2y d ^0( e ^-l) (A n (p) + A 22 (p)) 

= -2 / (e** - 1) £ - 2 • (32) 

J { ' «=i (po + \V ■ U) +(^ 2 + A« ) (p,c7,6)) 



The explicit integration over p yields 



. f <Pp eV'-l ^ F(p) 

The expressions of the eigenvalues A"^ (a = 1, 2) are given in Appendix B. 

In remains to specify the geometrical structure of the flow. For the sake of simplicity, we assume that the vector r 
corresponds to the polar axis and that the vector b is orthogonal to both r and U. It turns out that S 2 (r) can be 
expressed as the sum of two terms: the first one is associated with the null eigenvalue A^, while the second one 
depends on the nonzero eigenvalue A^ . Namely, 

S 2 (r) = -±(h(r)+I 2 (r)). (34) 



The expression of h(r) is derived in Appendix C: 



W-W/^^-Di^il (35, 



By simple algebraic manipulations, it can be recasted into the form 



Dn ^ T ( s+3+2n \ 2n 
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We can conclude that, for short distance r, the leading contribution in h(r) is r 2 , that is a dissipative contribution. 
Some lengthy algebra (see Appendix C for details) provides also an expression for ^(r)- 



E ■ 

i=l,2 



Em=0 s lmF m (x, p 2 ' S, S J + 



l3 



(l-a; 2 )3 



13 



IL* (ELo - *i0 ^ n^V; s, s)) 



TZ 2 J 



(37) 

where the coefficients and the functions F, are specified in Appendix C. 

The key remark for proceeding in this calculation is that the stochastic measure p 2 + s e~( Lp } dp gives a significant 
contribution to the first integral in (37) only in a narrow region of wavenumbers close to p where the function 
p2+s e -{L P ) k as its maximum, i.e. 



P= L 



1 s + 2 



Notice that the function 



jj 2 - p thus contributes to the integral by taking values close to K 
the sufficient condition (22) for the stability of small perturbations determines the upper bound 

TZ < 4(s + 2) , 



(38) 

4 ( a+2 ) Moreover, for p = p 



(39) 



This implies that for sufficiently small Reynolds' numbers the wavenumbcr p is stable. Under this condition, the 
leading contribution in (37) can be obtained by performing an expansion in powers of 8ly ^ p2 ■ 
One finally obtains the complete expression of the structure function (see Appendix C for details) 

1 



S 2 (r)=--(/i(r)+/ 2 (r)) 



+0 




n+ll 



s + 2n + 3 



1 + 5 



2^S 2n + 4 



r(2n + 4) S§ T(2n + 6) 



(i) ! 



for 1 < K < 4(2 + s) , 



(40) 



,4(2 + 5), 

At leading order in the distance r it is dominated by a dissipative contribution. 

We conjecture that this analysis can be extended to the parameter region defined by the conditon TZ > 4(2 + s), 
where the statistically relevant wavenumbers can be unstable. As shown in Appendix C, in this case gives 
two contributions: one is again dissipative, while there is another one yielding the nontrivial scaling behavior r 2 / 3 . 
Specifically, the expression of ^(r) for TZ > 4(2 + s) is 



Si{r)= 

TTV 



! + f 2 



n+l 



n=0 



r(§) W 



p ^s+2ri+3 N 

T(2n + 4) 
/3s + 3n + 5 



n=0 



(IT 



+ 



4(2 + s)' 



(41) 



This is dominated by the term r 2//3 for sufficiently small distances. Indeed, the crossover scale between the r 2 and 
the r3 terms occurs at 



r 
L 



FTZ- 



(42) 



In appendix C we evaluate the constant F ~ 0.6 and we report also the expression of the numerical coefficient Cb(E) 
(the general expression of the coefficients C„(S) appearing in (41), has been skipped, because it has no practical 
interest). 

It is a remarkable fact that S2 can exhibit the scaling behavior predicted by K41 theory, which is assumed to hold 
when the velocity fluctuations are turbulent in the so-called inertial range of scales. This suggests that hydrodynamic 
fluctuations in a system at the very initial stage of instability development already contain some properties attributed 
to the developed turbulence regime. 



7 



VI. CONCLUSIONS 



In this paper we have exploited field-theoretic calculations to reformulate the random forced Navier-Stokes problem 
in terms of a quadratic action functional. At a formal level, the latter has the same structure of the action describing 
thermal fluctuations in irreversible stationary processes. The crucial step for obtaining the hydrodynamic evolution 
operator which appears in the action functional, is the integration over all longitudinal components of both velocity 
and associated auxiliary fields. With respect to the standard formulation which is the starting point for diagramatic 
strategies, we thus perform one more field integration. The positive definite kernel in the action functional appears 
in the form of the inverse of the forcing correlation function. 

In terms of the action functional, the knowledge of the whole velocity statistics reduces to the computation of 
functional integrals. However, finding an explicit solution for these integrals is quite a difficult task, due to the 
nonlinear character of the problem. This forces us to introduce some approximations. The starting point is the 
identification of a stationary solution around which we linearize the evolution operator. We define also a fluctuation 
field (with respect to the stationary solution) and we are able to compute (perturbatively, in the inverse of the 
Reynolds number) the functional integrals over such fluctuation field. 

In principle, the strategy might be applied to evaluate any velocity multipoint statistical quantity. In order to 
reduce the complexity of the algebraic manipulations we limited ourselves to the calculation of the two-point second 
order momentum of velocity. Indeed, we aim at understanding if fluctuations at the early stage of their development 
(accordingly, we dub them as pre-turbulent fluctuations) already contain some important features of developed tur- 
bulence. We are interested, in particular, to scale invariance. In this respect, we find that fluctuations are organized 
at different scales in a self-similar way. Remarkably, the scaling exponent coincides with the dimensional prediction 
of the Kolmogorov 1941 theory [6] valid for developed turbulence regimes. Whether or not such exponent is a genuine 
reminescence of the developed turbulence phenomenology needs further investigations. 

Unfortunately, the complexity of the derivation leading to the K41 scaling law does not allow to identify precisely 
the very origin of such a dimensional prediction. We can however argue a relationship between the observed dimen- 
sional scaling and the conservation laws (for momentum and energy) associated to the two eigenvalues of the matrix 
appearing in the action functional (25). 

Another point to be emphasized is that the pressure term does not play any role in the derivation of the dimensional 
scaling law. This is just a consequence of the fact that all the longitudinal degrees of freedom can be averaged out 
from the very beginning of the computation. 

We want to conclude by outlining some open problems and perspectives. A first question concerns the relevance 
to be attributed to the solution around which we linerize the evolution operator. On one side we do not see any 
rigorous mathematical motivation for invoking the need of a unique solution. Just heuristic arguments based on 
physical considerations allowed us to identify the selected solution. Indeed, it represents a shear-like solution, which 
is a well-known generator of instability towards smaller and smaller scales. 

Another interesting point concerns the computation of the third-order moment of the velocity correlators. In this case 
the predictions of our approach could be compared with the 4/5-law, which is one among the very few exact results 
of turbulence theories. 

Finally, the extension of our results to other classes of transport problems, including passive scalar advection, could 
provide a better understanding of the basic mechanism at the origin of the observed scaling behaviors. 
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APPENDIX A 



In this Appendix we study the stability analysis of the solution V a relative to the linearized equation (20). We 
have already observed that V a is a quasi-steady solution for a time i<rc= . The Fourier transform of eq.(20) 
with respect to the x spatial variable yields: 



vk 2 u a 



t,k) + 



-k-Uu a 
2 



t,k) + 



1 

— e t d 
4 



U^k, 



t,k-C 



bAU 
AvTl 
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-u a It, k + C 



bAU 
AvTl 



+ U a C 



AvTZ 



AuTZ I \ ' AuTZ 
= . (Al) 
By performing a perturbative expansion up to second order in the parameter 7Z _1 , one obtains the system of equations 
d 



-uf 0) (t, k) + ^ 2 «f 0) (t, k) + l -k- Uuf 0) (t, fc) = 



^uf 1} (t, k) + uk 2 u a (l) (t, k) + l -k- Uufa (t, fc) 



-k-UC 
2 



4i/ft dk fj (°) 



i™( ?a£ 0*- 



4i/fc (°) 



(A2) 



(A3) 



2 fc ^ G AvTZ dk U 



I- ■ Uuf 2) (t, k 



1 (bAU) 
u h (t,k)--U-C^-^u^(t,k) , 



(A4) 



This system of equations yields the perturbative solution 

fi" (t, fc) = e-(" fca+ ^- £ )*(i^ 0) (fc) + Ffr (fc) 



k ■ U 



8i/£ 

- (Sac/) ■ jfej^j p) ^ 2 



J7 Q 

ft- c7 



6 A *7 ) • F (0) ( k ) t 



+ 



K 2 



(A5) 



For the expansion (A5) to be meaningful, the time t must be smaller than <~ 1Z. This amounts to impose the condition: 

8vK 2t d 



0<t< 



u 2 n ' 



(A6) 



This stability condition implies that, at any time t, the linear term in the curly brackets cannot overtake the exponential 
factor. Such a requirement can be traslated into the following spectral condition 



8v 2 K 
~U 2 ~ 



k 2 >l. 



(A7) 



APPENDIX B 

In this Appendix we sketch the calculation of the eigenvalues of the matrix M^(p). We exploit a perturbative 

approach, whose expansion parameter is ^. The matrix M^(p) acts on the two-dimensional space of the transverse 
functions and on the one-dimensional space of the longitudinal functions. Only the transverse degrees of freedom are 
physically meaningful. 
In terms of the ^ expansion we have 

M = M {0) + M (1) + ... (Bl) 

where 
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M (o) 



i[Po + ^P-U ) + vp 2 



C 



A complete orthonormal basis in R 3 is given by the vectors 



(bAuV r (^ At7 ) 



4vK 



(B2) 



irf= 



6 Ap 



g(p) (bApY -f(p) (UAp)' 



^f7¥Wf(p)Kp)-9 2 {p) ' 



n?= 



p" 
p 



where we have defined 



f{p) = b 2 p 2 - (b ■ p)\ g(p) = (b ■ U)p A -(b-p)(U-p), h{p) = U z p z -(U-p) 



(B3) 



(B4) 



Here, II" and n 2 span the transverse subspace, while II3 spans the longitudinal one. Likewise, the eigenvalues can 
be represented in terms of a perturbative expansion as 



X a = \ a {0) + X a {1) + ... where a =1,2, 3. 
The zero-order eigenvalues A" ^ are degenerate and have the form 



A (0) 



* ( Po + ^P ■ U ) + vp 2 



(B5) 



(B6) 



The evaluation of the first order corrections X?^ requires the diagonalization of the matrix with elements = 
(Ilj, M^IIj), — 1,2,3). After some simple but lengthy calculations we find 



'(1)21-^(1)12 



*(i)= \ y M (i)ii + M(i )2 2 - y / (M (1)11 +M (1)22 ) 2 +4Af ( 
A(i)= \ (m (1)11 + M (1)22 + yj (M (1)11 + M (1)22 ) 2 + 4M (1)21 M (1)12 ^ , 



A^)= M(i) 33 



(B7) 



and 



M (1)11 = 



C 



■P 



-w(p) , 



Mm 129= — 



16uH f(p) 

c ((bAu)-p)(b-p)g(p) 



M (1)12 = 
M (1)21 = 

M(i) 33 = 

where we have defined 



(1)22= ~16^ f(p)(f(p)h(p)-g 2 ( P )) 
C ((bAij).p)(b.p) 
1 ^.f(p)Vf(p)Hp)-9 2 (p) 
C ((bAU)-p) 

1 ^f(p)VfiP)Hp)-9 2 (p) 



(p-u)w(p) + (b-U)g(p)-U 2 f(p) 
(b- U) g(p) + 2 (p- U) w(p) - U 2 .f(p) 
(b-u) [(b-p)g(p)- (p-U) f(p) 



(B8) 
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w( P )=b 2 (p-U)-(b-p)(b-U) . (B9) 
With the particular choice performed in Section V, the two physically relevant eigenvalues are 

Aji)=0, 

A(i)= Yoj/ji { sm ^ cosOu [cos 2 4>u + cos (2(<f>u - 4>))] sin 2 9 

+ cos 2 9 V sin 29 cos(<f>u - 0)} . (BIO) 



APPENDIX C 

This Appendix contains the essential steps necessary for computing the second order structure functions. 
We start from the expression of S 2 



S(r)=- f ^ e * P ' P ~ 1 V F{P) 



= -\{h(r)+h{r)) , (CI) 

By considering the explicit expressions of the statistical function F{p) and of the eigenvalues A"^ (see eq.(BlO) ), one 
has 



AM-W/^e"-') 1 *^. (02) 



d 3 p (e r P- p - 1) {Lp) s e-^ 

r) 3 p 2 + 16 ^ TC [2 sin 9\j cos #;y sin 2 cos 2 (p + cos 2 #;y sin 29 cos < 



(03) 



By exploiting translational invariance, in (C3) we have applied the transformation 4>u — <p ^ —<j)- The evaluation of 
(C2) follows from a standard procedure: 



2tt 2 £j (2n)!(2n+l) 7 ^ ^ 

= _^ r 2 y (_i)«+i Li^^lfM 2 " . 



(2^) 2 ^ ' T(2n + A) \L 
Concerning the term hfr) we first perform the integration in the <j) variable. Namely, 



(C4) 



I 2 (r) = D L 3 J™ 0^(Lp) s e- (L rf J ^ d(coa 9) (e^ rcos9 - l) J (C5) 



where 



Jo 



p 2 + 16 ^ TC [2 sin cos 9jj sin 2 6* cos 2 <\> + cos 2 #[/ sin 29 cos 0] 

32i/ 2 K / zdz 

U 2 L az 4 + bz 3 + cz 2 + bz + a ' ( ' 



= —i 

i - 

'7 



with z = e 1 ^ and the integration is on the unitary circle. The coefficients a, b, c are given by 



11 



a = sin 9u cos 9\j sin 2 9 , 6 = 2 cos 2 9u sin cos 9 

2,2v 2r R 2 „ . „ „ . 2 „ 
c = — —5 — p + 2 sin 0[/ cos 9u sin # . 
t/ 2 



(C7) 



The evaluation of the integral (C6) requires the knowledge of the root of a fourth degree algebrical equation. By 
exploiting the Euler method [9] we end up with the expression 



Zi Zi 1 x, 



Sv 2 Tl 



p 2 ;£,~ 



X3 



(l-x 2 )* u=o 
The following definition has been adopted: 

8u 2 TZ 



%v 2 Tl 



X3 



(l-x 



2\3 



i = 1,2,3,4 . 



Fi=Fi \x, 
"I 12 



U 2 



A 



4 (l-x 2 ) 
64 1 - x 2 8v 2 K 



81 v2 x 
2 2 (l-x 2 ) 



90 



£ 2 x 2 



?7 2 



r 189— 



382 



120— 



E (l- x 2) 2 S(l-x 2 ) E x 



U 2 



-P 



504 



47 S 2 



S 2 (l-x 2 ) 2 E 2 x 2 (l-x 2 ) 



%v 2 n 



u 2 



-p 



32 S 2 



S 3 x 2 (1 - x 2 y 



x [<• 1 + (1 - 4 x 27 h)* 
1 (l-x 2 ) 1 %v 2 K 2 



1 - (1 - 4 x 27 h)- 



4 4 

1 £3X3 



(l-x 2 ) 3 



3 S 



3X3 



u 2 



3£§S xl (l-x 2 )? 



with 



X = COS , S = sin #[/ COS 6*;y , S = COt 6*(7 

/111 




Here e is the cubic root of the unity: e = Mjr^ - Finally the function ft, is 

^v 2 Kp 2 17£ 2 x 2 + 3(l-x 2 ) 



16 + 30 £ J 



x 



111_ E 4 ^ 



16- 



+48 



8v 2 1Z ,r • 2 



t/ 2 h ; (1 



l-x 2 4 (l-x 2 ) 2 

3 

x -128 + 81 £ 4 



1 



U 2 E 

x 4 



(l-x 2 ) 2 
81. 



(l-x 2 ) 2 



8v 2 TZ p 2 



-180E Z 



+ 



l-x 2 U 2 E\ (l-x 2 ) 3 



2^,2 



378- 



+ 764 



S 2 x 



+ 



8^ P 2X 2 
U 2 E 



Y?x 2 

1008- — + 94- 



(1-x 2 ) 3 (l-x 2 ) 2 



(l-x 2 ) 2 
+ 64 



240- 



(1 



„,2^3 



l-x : 



8v 2 K P 2s 3 



U 2 E) (l-x 



2V3 



Only the roots z\ and Z2 are included into the unit circle, therefore (C5) becomes 



(C8) 

(C9) 
(C10) 



(Cll) 
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I 2 {r)=D L 

*E 



Em=0 s imF m (x, ^jjf^p 2 ; E, ^ + iS—^ 



f7 2 

(I-- 2 ) 



(l-iE 2 )? 



= 1,2 



13 



(ELo {sik - s ik ) F k (x, ^pi; S, S)) 



(C12) 



As we have already observed in Section V, only the values of the variable p around P = j;y ^-^r give a significant 
contribution to the integral in (C12). We observe that ^jjr^p 2 — ► 4 ^ 2 ^ and the stability condition (22) imposes: 

Kft < 4(s + 2) . (C13) 

The evaluation of the leading terms is then possible by performing an expansion in the parameter ^ n P~ 2 — * -f^ -2 
that, by virtue of (C13), is smaller than unity if ( < 

For C > \/ t ne contribution to the integral rapidly vanishes. For 1 < 1Z <C 4(s + 2) we obtain 



S 2 (r)=--(/!(r)+/ 2 (r)) 



D 



(2tt) 2 ^ 



E (-!)" +lr 



n=0 



s + 2n + 3 



1 + S 2irS 2n + 4 



r(2n + 4) E§ r(2n + 6) 



2n 



+0 



ft 



4(s + 2) 



(C14) 



By extending the validity of our calculations to 1Z > 4(s + 2), we have ^^p 2 — > t|-C 2 < 1 for £ < y ^t 2 -. As in the 
previous case, we expand (C12) in power of the parameter ^-( 2 < 1 and we obtain: 



J 2 (r)~ SL> 0j L 2 



+^E 



(27T) 2 

(l-z 2 ^ 



C'e - ^ / dx (e l ^ x - 1) 



1 + ^ + 



E m =o (ar,0;E,S) 



8^2 9 



Em=0 ( X > 2/! S : S ) 



(C15) 



Two different terms, / 2 (r) = I^i 7 ") + -^f( r )> can be identified in (C15). The evaluation of the first is straightforward 
and we obtain 



5£> 2^ (-l) n+1 



r 2 E 



2(2^) 2 ^r(2n + 4) 



s + 2 + 2n 



7^ 



s + 5 + 2n\ \ / r 



(z)' 



x 1 + 



1 4( S + 2) 



ft 



The evaluation of the second term is more cumbersome. The leading term can be recasted in the form: 



B 8DpL 2 [°° d( _ C 2 



1 — x 



2\ 3 °° 



X3 



E^(^ 2n 



(C16) 



(C17) 



n=0 



The coefficients A t are E-dependent numerical constants. The first two of them are given by the expressions 
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4>(E)= 



16\/3 (l - sin |) cos (| tan" 1 V26) 



^i(S)=-- 



65sin(|tan- 1 % /26) 



512V26cos 2 (itan" 1 ^) 

The exact form of these coefhcients is however irrelevant for our analysis. After some calculations we obtain 



3s + 3n + 5 
6 



(z)' 



*r(!) 

where the coefficients C„(S) are depend on the constants A t . For n = we have 

v 7 27V3 9V3 
Comparing the scaling behavior of the term /.f (r) with the term /^(r) we find 



2 x 8.328-v/tt- 



0.0336- 0.1127 cot 2 T 



2 + sin 0(7 cos 9u 



(C18) 



(C19) 



(C20) 



(C21) 



For the expansion in term of to be meaningful, the parameter 9u must lie in a small region around the value \. 
This implies: 



r ~ FTT^L, with F ~ 0.6 . 



[1] P. C. Martin, E. D. Siggia, H. A. Rose, Statistical Dynamics of Classical System, Phys. Rev. A8 (1973) 423 

[2] C. Dc Dominicis and PC. Martin, Phys. Rev. A 19, 419 (1979); J. P. Fournier and U. Frisch, Phys. Rev. A 17, 747 (1978); 

V. Yakhot and S. A. Orszag, J. of Sci. Comp. 1, 3 (1986). 
[3] L. Ts. Adzhemyan, N. V. Antonov, and A. N. Vasilev, The Field Theoretic Renormalization Group in Fully Developed 

Turbulence (Gordon & Breach, London, 1999). 
[4] L. Bertini, A. De Sole, D. Gabrielli, G. Jona-Lasinio, C. Landim, Fluctuations in Stationary Nonequilibrium States of 

Irreversible Processes, Phys. Rev. Lett. 87 (2001) 1 
[5] L. Onsager and S. Machlup, Phys. Rev. 91, 1505 (1953). 

[6] A. N. Kolmogorov, Dolk. Akad. Nauk SSSR 30 (1941) 9; Dolk. Akad. Nauk SSSR 31 (1941) 538; Dolk. Akad. Nauk SSSR 
32 (1941) 16; U. Frisch, Turbulence; the legacy of A. N. Kolmogorov (Cambridge U. press 1996) 

[7] L.Ts. Adzhemyan, N.V. Antonov and A.N. Vasilev, Renormalization Group, Operator Product Expansion, and Anomalous 
Scaling in a Model of Advected Passive Scalar, Phys. Rev. E 58 (1998) 1823 

[8] C. Kipnis, C. Landim, Scaling Limits of Interacting Particle Systems (Springer, New York, 1999). 

[9] W. S. Burnside and A. W. Panton, The Theory of Equations - with an introduction to the theory of binary algebraic forms 
- V.l (Dover Publications, inc. New York 1912) 



14 



